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Abstract 

Second order circularity, also called properness, for complex random variables is a well known and studied concept. In the 
case of quaternion random variables, some extensions have been proposed, leading to applications in quaternion signal processing 
(detection, filtering, estimation). Just like in the complex case, circularity for a quaternion-valued random variable is related to the 
symmetries of its probability density function. As a consequence, properness of quaternion random variables should be defined 
with respect to the most general isometries in 4 D, i.e. rotations from 50(4). Based on this idea, we propose a new definition 
of properness, namely the (p, 1 , /r 2 )-properness, for quaternion random variables using invariance property under the action of 
the rotation group 50(4). This new definition generalizes previously introduced properness concepts for quaternion random 
variables. A second order study is conducted and symmetry properties of the covariance matrix of (p 1 , fi 2 )-proper quaternion 
random variables are presented. Comparisons with previous definitions are given and simulations illustrate in a geometric manner 
the newly introduced concept. 


I. Introduction 

Quaternion random variables and vectors have attracted a substantial amount of attention over the last decade in the signal 
processing community. Amongst the interesting problems arising in random quaternion signal processing, the extension of 
the concept of properness (also called second order circularity) of quaternion random variables has been investigated by 
several authors. Orignal work was made by Vakhania [1] who studied infinite dimensional proper quaternion vectors in right 
quaternionic Hilbert spaces. His work was transposed to quaternion random variables in [2], Several papers exploring the 
concept of properness for quaternion valued random vectors and signals, and their applications, followed soon after [3], [4]. 
Tests to check the properness of quaternion random vectors were proposed in [5], [6]. Quaternion properness has then found 
applications among which widely-linear processing [7], [8], detection [9], [10], ICA [11], [12], adaptive filtering [8], [13], 
study of random monogenic signal [14], quaternion VAR processes [15], Gaussian graphical models [16] and directionnality 
in random fields [17]. 

Quaternion signal processing makes use of the definition of properness given in [2], [3] which relies on the vanishing 
of some covariance coefficients between quaternion components. Properness can also be stated in terms of invariance of the 
probability density function under specific geometric transformations: Clifford translations ([18]). Indeed, left and right Clifford 
translations have been used in [3] and [2] respectively to define quaternionic properness. However, Clifford translations are not 
the most general isometric transformations in 4D space. They are actually very specific isometries. In this paper, we propose 
a new definition of properness for quaternion random variable, based on the most general isometric transformation in R 4 , 
namely the action of 50(4). This approach relies on geometric properties rather than arguments based on the vanishing or 
not of correlation coefficients between the components of the quaternion random variable. In particular, it is shown that the 
most general case of properness does not imply any vanishing components in the covariance matrix of a quaternion random 
variable, eventhought the probability density function of this variable already exhibits strong symmetry. 

The proposed definition of properness, named (/r ; ,/x ? ,)-properness, leads to extra classes of quaternion random variables 
that were not identified previously. Already known classes are special cases of the introduced definition. A second order study 
(covariance matrix analysis) illustrates this in the well known case of Gaussian quaternion random variables. It is demonstrated 
in this paper that (fi h fi r )-properness is the most general definition to study the geometry of quaternion random variables. 
An interesting property of the proposed approach is that, as it is not solely based on properties of the covariance matrix, it 
allows higher order symmetries study and open the door to potentialy new development in higher order statistics for quaternion 
signal processing. Actually, the (/r z , fi r )-properness study proposed here could be extended later to (/r ; , fi r )-circularity for 
quaternion random variables. 

The paper is organised as follows. In Section |IIJ we introduce some elements of quaternion algebra and results about 
the rotation group 50(4). Then, in Section [lH] we present our definition of properness for quaternion random variables. 
The properties of the covariance matrix of proper quaternion random variables are examined in Section |IV[ together with 

presents the conclusions and perspectives of the presented work. 

II. Quaternions and 4D Geometry 

After a rapid introduction to the algebra of quaternions, we present their close relation to the rotation group 50(4). 


comparison to existing work. Section M provides some examples of proper quaternion random variables and Section VI 
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A. Quaternions 

Quaternions are 4D hypercomplex numbers discovered by Sir W.R. Hamilton in 1843 [19]. The set of all quaternions is 
denoted EL A quaternion g £ HI can be expressed in its Cartesian form as: 

q = a + bi + cj + dk ( 1 ) 

where a, b, c, d £ ffi. are called its components. The three imaginary numbers i. j. k fullfil the well known relations: 

i 2 = j 2 = k 2 = ijk = -1 ( 2 ) 

The scalar part of q is S(g) = a and its vector part is V(g) = q — a. The three imaginary parts of q are also denoted 
Ry ( q ) = b , $ry (g) = c and (g) = d. The restriction of a quaternion q to a degenerate quaternion is itself a quaternion 
that is made of some components of q. For example, the restriction q[ij] = bi + cj = (q) i + Ay (q) j is only a 2D 

hypercomplex part of q made of its Ai and Ry components. Restrictions can be made using three indexes as well, with for 
example, q[ij,k] = bi + cj + dk a 3D hypercomplex sub-part of q. A restriction with one index is simply g^j = bi for example. 
The restriction gm = a is the real part of g, while gri,*,j,fe] is <7 itself. By extension, one can also make use the restriction 
notation for HI by noting for example that qpj.k] € H^yy.], with H^yy.] = M © j'K. © fcR. Subfields of HI isomorphic to the 
complex field will be denoted C tJ when = R CD pM. and p 2 = — 1. 

The product of two quaternions p, q £ HI is not commutative, i.e. qp 7 ^ pq; see for example [20] for the expression of the 
product components. The conjugate of a quaternion q is denoted q* and defined as: q* - a - bi — cj — dk. The modulus 
of q £ HI is denoted \q\ and given by: |g| = (a 2 + b 2 + c 2 + d 2 ) 1 ^ 2 , while its inverse is: g -1 = g*/|g| 2 . A quaternion g is 
called pure when its scalar part is zero, i.e. S(g) = 0. The set of pure quaternions is denoted V(HI). Quaternions with modulus 
equal to 1 are called unit. Unit quaternions form a group under the quaternion multiplication. This group is sometimes denoted 
Sp{ 1), and can be identified with the unit sphere S 3 C R 4 . In HI, in addition to conjugation, it is possible to define involutions 
[ 20 ] the following way: 

Definition 1: Given a quaternion q € HI and a pure unit quaternion p £ V(H), then the mapping q —> g" M such that: 

q‘ u = -MU ( 3 ) 


is an involution. 

Properties of involutions that will be of use later in the paper are: 


pq IJ = p^q^ 

(4) 

q-n 11 = q-(nn) 

k rr = 

for q,p £ HI and p,ij £ V(HI) such that {1, p, rj, prf\ is a basi^jin HI. Note that g" M consists in a rotation of V(g) of an 
angle 7 r around p. Any quaternion q £ HI can be expressed in its Euler form as: 

q = \q\ (cos 0 q + p q sin 0 g ) = |g|e AX « 6 ' 3 (5) 


where: 


0 q = arctan 


^ \/b 2 + c 2 + 


bi + cj + dk 

^ 9 Va 2 + b 2 + c 2 + d 2 

p q is called the axis of g and 9 q is called the angle of g. If g is a unit quaternion, then it can simply be expressed 
The set of unit quaternions consists in the 3-sphere in R' 1 : 


( 6 ) 

(7) 

as e fx i eq . 


S 3 = {geH||g| = 1} ( 8 ) 

A pure unit quaternion takes the form c; , R 2 = p q . Also, any pure unit quaternion p has the property that p 2 = —1. An other 
important way of looking at quaternions is to consider them as pairs of complex numbers. This can be nicely done using the 
Cayley-Dickson form of a quaternion q = a + bi + cj + dk which reads: 


q = Zi + z 2 j 


'The set { 1 , fi, 77, p .77} is a basis in HI if p 2 = 77“ = 


1 and /x_l_?7 = 0, i . e . /X77} = 0. 
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where Zi,Z2 £ C are given by z\ = a + bi and Z2 = c + di. This notation will be of use when studying quaternion 
random variables in Section [III] In addition to Euler and Cayley-Dickson notations, we mention the matrix representation of 
quaternion product as it will be used later on. First, recal that as a vector space over R, H is isomorphic to R 4 . Any quaternion 
q = a + bi + cj + dk £ H can be represented by the real vector qm = [a,b,c,d] T . Then, the quaternion product qp of 
q = a + bi + cj + dk and p = po + Pi* +P2J + P 3 ^, denoted L q {p) = qp can be represented by the matrix-vector product: 


L g(p) 


a —b —c —d 


Po 

b a —d c 


Pi 

c d a —b 


P2 

d —c b a 


P3 


(9) 


This is a left quaternion multiplication of p by q, which explains the notation L q {p). Now, for a right quaternion product of 
p by q, i.e. R q {p) = pq, one gets the following matrix multiplication expression: 


R q {p) 


a —b — c —d 


Po 

b a d —c 


Pi 

c — d a b 


P-2 

d c —b a 


P3 


( 10 ) 


As noted in [21], the sets of matrices L q and R q as debited above, are two subalgebras of Ml (4. R j^] both isomorphic to EL 
An interesting property of the matrix multiplication by L and R is as follows. 

Property 1: Given any two quaternions q, r £ H, then: 


L r {Rq(p)) = Rq (L r (p)) = rpq (11) 

for any p £ EL 

This means that the matrix representations of right and left products commute. In addition to this property, note that any matrix 
M £ A4(4,R) can be expressed as a linear combination of R a Lf, matrices [21]. 


B. Rotations in 4D 

It is well known that quaternions are efficient at representing 3D and 4D rotations [18], [22]. This ability of quaternions 
to encode rotations is due to their relation to the Spin groups Spin( 3) and Spin( 4) [21]. We give here a rapid overview of 
the theory of rotations in R 4 using quaternions. Material about rotations in 4D space can be found in various forms (matrices, 
quaternions, Clifford algebras) in [22], [21], [18], [20], [23]. First, we recall some known results about 50(4). 

1) The rotation group 50(4): The group of rotations in R 4 , denoted 50(4), is a matrix Lie group. An element of 50(4) 
can thus be represented by an element of A4(4,R) . 

Definition 2: A matrix U £ A4(4, R) is a rotation matrix, i.e. an element of 50(4), iff: 

UU T = I 

det(O) = 1 

where I is the 4x4 identity matrix. 

As 50(4) is a matrix Lie group, any matrix U £ 50(4) can be expressed as [24]: 

U = e s (14) 

where E £ A1(4,R) is anti-symmetric, i.e. E T = —E. As E is anti-symmetric, its eigenvalues are purely imaginary and come 
in pairs, say ±.ia and ±i/3. As a consequence, the eigenvalues of the rotation matrix U £ 50(4) are e~ za , e _ * /3 }. 

It is well known that rotations in R 4 are either simple or double rotations [23]. A simple rotation leaves an entire 2D plane 
invariant, while a double rotation leaves only the center of rotation O invariant. Simple rotations can be of two types: left 
isoclinic or right isoclinic. It is also a well known fact that the set of left isoclinic rotations form a non-commutative subgroup 
of 50(4) that is isomorphic to the multiplicative group ot unit quaternions S 3 . In a similar way, the set of right isoclinic 
rotations is isomorphic to S 3 . A simple rotation is parametrized by an angle and a 4D vector (which debnes the plane left 
invariant by the rotation). A double rotation in its general form is given as follows. 

Property 2: A double rotation matrix U £ 50(4) can be decomposed into the matrix product: 

U = L a Rp (15) 


( 12 ) 

03 ) 


2 A4(4,R) is the set of 4 X 4 real valued matrices. 
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where L a is a left-isoclinic rotation with angle a and Rp is a right-isoclinic rotation with angle 8. 

Note also that, thanks to Property [I] U = L a Rp = RpL a as left and right isoclinic rotation matrix representations commute. 

The two angles of the double rotation are simply the eigenvalues of U mentioned previously [21], The action of left and 
right isoclinic rotations are illustrated in Figure [T] displaying the action of these rotations in 4D by their simultaneous actions 
in two non-intersecting (completely orthogonal) 2D planes in R 4 . 




Fig. 1. Rotations induced by a left isoclinic rotation (two upper planes) with angle a and by a right isoclinic rotation (two lower planes) with angle /3 in 
two completely orthogonal 2D planes of R 4 . 


A double rotation is thus defined by two axes (4D vectors, to be specified later), and two angles a and /3 (the eigenvalues of 
the associated rotation matrix). The cases where the ratio between the two rotation angles is irrational lead to rotations e* >: f I 
for all t > 0 and t £ R [21], Only completely symmetrical objects in 4D (4D hyperballs) are left invariant by such rotations 
with irrational value of /3/a. In the sequel, we will only consider cases where /3/a = k with k £ Q, leading to rotations in 
4D that describe the symmetries of objects invariant under discrete 4D rotations. 


C. The group Spin{A) 

Pairs of unit quaternions (u, v) give the product group 5 3 x S 3 . This group acts on H ~ K 4 as: 


Q 


uqv 


(16) 


where q £ H, and as it preserves length and is linear, each element (u. v) corresponds to en element of SO (1). One can also 
see from eq. © that (—it, —v) corresponds to the same element of 50(4), and from this it comes that 50(4) is the group 
S 3 x S 3 with ( u , v ) and (— u, —v) identified. This reads: 


50(4) 


5 3 x S 3 

{( 1 , 1 ),(- 1 ,- 1 )} 


(17) 


where S 3 x S 3 = Spin{ 4) is the Spin group: the Lie group that double covers 50(4). 

From eq. ( |~i~ 6 | ) and © it can be seen that the 4D rotation expressed with quaternions is parametrized by two unit quaternions. 
Given a quaternion q, left multiplication by a unit quaternion u consists in a left isoclinic rotation of angle 9 U and right 
multiplication by a unit quaternion v consists in a right isoclinic rotation of angle 9 V . In order to know in which planes 
are occuring the rotations, one needs to refer to the axes of the unit quaternions. As in eq. © consider performing a left 
isoclinic rotation of quaternion <7 by u = e M " a " and a right isoclinic rotation by v = First, decompose the quaternion 

q = a + bi + cj + dk using two of its possible Cayley-Dickson form^J 


q =z 1 + z 2 v u = (o' + b'n u ) + (c 7 + d'p, u )v u (18) 

q =wi + w 2 u v = ( a" + b" n v ) + (c" + d!' \i v )v v (19) 


where {1 , p, u ,i/ u , fi u i/ u } and {1 , /j, v ,v v , fj, v is v } are bases in EL See [20] for more details on Cayley-Dickson forms of 
quaternions. Note also that for any Cayley-Dickson form such as those in eq. © and ©, zi and z 2 (respectively w\ and 
w 2 ) define two completely orthogonal planes in R 4 (only intersecting at the origin). The left isoclinic rotation by u will perform 
a 9 U rotation in the zi -plane, together with a 9 U rotation in the 22 -plane. A similar reasoning with right isoclinic rotation as 
a right multiplication of q but the unit quaternion v = e* 1 " 6 ’’ leads to the conclusion that a 9 V rotation is performed in the 
u'i -plane while a rotation of —9 V is performed in the W 2 -plane. These rotations in orthogonal planes in 4D are displayed in 

Fig- [[] 


III. (/!;, /X r ) -PROPER QUATERNION RANDOM VARIABLES 

The extension of the concept of properness from complex to quaternion random variables is possible through the gener¬ 
alization of the notion of circularity for 4D random variables. From a geometric point of view, circular quaternion random 
variables have probability distributions exhibiting symmetries in R 4 . Before focusing on properness, we introduce the more 
general concept of circularity over EL 

3 Given a quaternion q G H. there is an infinite number of Cayley-Dickson forms, associated to all the possible basis of H. 
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A. Circularity 

In the case of complex random variables, it is possible to give a definition of circularity by means of characteristic function 
or moments symmetry properties for example [25]. An equivalent, and may be more intuitive and geometrical way is to use 
invariance of the probability distribution [25]. Here, we make use of this last mean to extend circularity to the case of quaternion 
random variables. 

Definition 3: A quaternion valued random variable q £ H is called circular iff: 


q = Iqr 
= e^qe^P 


( 20 ) 

( 21 ) 


for any pair of unit quaternions l = e Mi “ and r = i.e. for any l,r £ H and |Z| = |r| = 1 with respective axes fi l and fi r 
and respective angles a and 8. Notation = stands for equality in distribution. 

Definition |3] is very general as the angles of a and 8 both span [0,27r] and axes fi l and /x r are elements of S 2 without any 
restriction. Notation used in (20 1 is classicaly used to name symmetry groups and symmetrical objects (polytopes) in R 4 . See 
for example the naming of 4D chiral groups by Conway in [22], As mentioned previously. Definition [3] is geometrical , and 
means that the probability distribution of q is invariant by any double rotation from SO(4) made of the left and right isoclinic 
rotations represented by l and r. A circular quaternion random variable thus possesses a distribution invariant by any rotation 
in 4D, i.e. a distribution completely symmetrical with respect to the origin. 

In defintion[3] axes //, and /j, r are allowed to take any value over S 1 . Forcing these two axis to belong to a common basis 
in H (and thus to be orthogonal to each other) allows to define a lower degree of circularity for quaternion random variables: 
the /x 2 ) / 3 -circularity. 

Definition 4: A quaternion valued random variable q is called °(p. t , [i f 1 -circular iff: 


d 

9 = 


e^ a q e ^P 


( 22 ) 


for any angles a and /3, and provided that {1, fi 1 , /x 2 , /x 3 } is a basis in H, i.e. D 1 D 2 = ~ H 2 H 1 = Ha- 

Recall that the axes used in Definition [4] should be simply chosen such that // , _L fi 2 in order to have a basis of H. In the 
sequel, we will make use of basis {1, /j, 1 , fi 2 ■ /tt 3 } to study properness. Note that it is always possible to express the presented 
results in the standard basis in H, namely {1, i,j , k } through a simple change of basis, i.e. a rotation in R' 1 . 

Note also that the first following restriction on angles a and 
a = k8[ 27 t | 


in eq. ( | 22 ] > is assumed: these two angles should be multiple 
this condition ensures that not only completely 
symmetrical 4D distributions (i.e. circular ones) can be described with the invariance given in eq. ({22]). The second restriction 


of each other. 


with k € Q. As explained in Section II-B 


is that we will only consider second order circularity in the sequel. This means that a and (3 will simply equal 0 or 7 t/ 2 . 
This is motivated by the fact that we will focuse on properness of Gaussian quaternion variables, for which only second order 
statistics, i.e. covariance, needs to be considered. 

A complete study of /x 2 ) /3 -circularity for other values of a and (3 is out of the scope of this article and is left for 


future work. In the sequel, we focuse on second order 
“(/x 1; /x 2 )^-properness. 


*(Hi, // 2 ) -circularity for quaternion Gaussian random variables. 


B. a {ni, n 2 )^-properness 

Previous definitions of quaternionic properness in the the literature have only used left isoclinic rotations [3] or right isoclinic 
rotations [2]. As mentioned previously, these are not the most general rotations in R 4 . We now introduce a new definition for 
properness of quaternion random variables. 

Definition 5: A Gaussian quaternion valued random variable q £ H is called a (p ll p 2 )^-proper iff: 


d 

q = 


e »i a qe »20 


(23) 


with a, (3 € {0, 7 t/ 2 , — 7 t/ 2 } and {1, /x l3 p 2 , HiH 3 } a basis of EL 

The fact that we focuse on particular angles allows a simplification in properness denomination/study. In the sequel, we will 
use the following simplifed notations: 

• (1,/x 2 )-proper: if a = 0 and fi = 7 t /2 

• l)-proper: if a = 7 r /2 and (3 = 0 

• {p 11 /r, 2 )-proper: if a = n/2 and j3 = 7 t /2 

This notation is unambiguous in the special case we are looking at, and it avoids superscripts. However, it should be kept 
in mind that the general notation is necessary in the more complicated study of circularity. We will thus use the (p 11 p, 2 )- 
properness notation to study Gaussian proper quaternion random variables in the sequel. 

The proposed Definition [5] of properness has some properties. We now give some of them which are to be used later. 
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Property 3: Given two pure unit quaternions // and v and a set of 2N angles ai,... ,ajv toid [3 \,... ,(3n, a quaternion 
Gaussian random variable which is, Vi, “‘(/r, v) 1 ' -proper, is also v)' -proper, with £ = YliLi a i and X = Xa=i Pi- 
Proof: By direct calculation for N = 2, using the fact that e M “e M0 = e ^ a + e ) when exponentials share the same axis. 
The proof follows easily by induction. ■ 

Property 4: Consider two angles a et /3 equal to ±7 t/ 2 and 2 N pure unit quaternions /tt 1 , /x 2 ,..., p N and u \ .v 2 ,.... r'v 
from a common quaternion basis. A quaternion random variable which is a (p i , i/q)^-proper Vi € 1,2..... A' is also a (rj, £,)'- 
proper with : 

N 4— 

V=Y[p-i= P-nPn-i ■ ■ • H2M1 

< i=1 
->N 

€=n Vi= Vih>2 ■ ■ ■ un ~ i v n 

< i= 1 

and where axis p and £ belong to a common basis of H. 

Proof: Using the fact that for a pure unit quaternion p, one has t ±IL7T ' 2 = ±p, and respecting the order in right/left 
multiplication, the proof is completed. ■ 

Property 5: If a quaternion Gaussian random variable q is ±7r/2 (/r. p) ±; ' 2 -proper with p £ {i,j,k}, then it is as well 
both ±lx ^(y, £) ±7r/, "-proper and ±7r ^ 2 (£, i/) ±7r,/ “-proper if is a quaternion basis. 

Proof: By simple calculation. ■ 

Just like properness in the complex case, (p 1 , p , 2 )-properness induces some structure in the covariance matrix of a quaternion 
random variable. We study in the next Section the consequences of (p 1 , /r 2 )-properness on the covariance matrix of a Gaussian 
quaternion random variable q. 


IV. Covariance matrix of (/r 1 , /t 2 )-proper Gaussian quaternion random variables 

It is well known [2] that in order to study the covariance of a quaternion random variable, one can choose amongst the 
three equivalent vector representations : 

• = [a, b, c, d] T , q R £ R 4 

• <7c = [zuz*, z 2 , z \] T , qc £ C 4 

• m = ®eD 4 

where * represents complex conjugatiorj^] and 7f u with v = p, t , p, 2 . p, :>t are the involutions introduced in eq. (|3j. The complex 
notation is assumed to be q = Z\ + z 2 p 2 with Zi,z 2 £ C /i . | . Note that the choice for these three vector notations is arbitrary as 
it depends on the basis chosen in H. Here, we make use of the basis {1, /x l5 p 2 , ;j, :( } in H to be consistent with the properness 
definition in eq. (23 1 . This will allow comparison with existing work [3] in Section |IV-C 


A. General definitions 

We begin with general definitions for quaternion Gaussian random variables, their covariance matrix and their natural 
symmetries. 

Definition 6: A centered quaternion valued random variable q is called Gaussian with covariance matrix (in its quaternion 
representation Qh) T q .H, i-e. Q ~ A r (0, k q .i-0, if 4s probability density function takes the form: 


pq(qh) = 


y/2n det(T q , H ) 


g 2 


(24) 


with the explicit expression of the covariance matrix T q H given in (25 i and where det (r qi ni) is the (/-determinant of T q H as 
defined in [26]. 

The covariance matrix of a quaternion Gaussian random variable is given using the following definition. 

Definition 7: Consider a centered quaternion random variable q, i.e. E[g] = 0. The covariance matrix of <y, using its quaternion 
vector notation qjj, is: 


r q , H = EfqnqJj] = 


(T 2 

7i/x, 

7in 2 

7im 3 

7mi1 

a 2 

7/x 1 /i 2 

7m ,/X 3 

7m 2 i 

7rt 2 Mi 

a 2 

7/X 2 (i 3 

7/Xa 1 


7#i 3 /x 2 

<7 2 


(25) 


Conjugation is not ambiguous here as we consider real quaternions and complex numbers. Conjugation thus only consists in negating the imaginary part 
which is 3-dimensional for quaternions and 1-dimensional for complex numbers. The notation * is thus used over H and C. 
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where f stands for conjugation-transposition and with the use of the following notation: 


E[M 2 ] 

= a 2 

E [?rn 

= 7i^ 

E [q^q*] 

= 7Ti 

Et^rn 

= Irju 


(26) 


and with r\,v £ V(H) being unit pure quaternions of the chosen basis in H, i.e. taking values p .,, /i, 2 and // 3 . 

As r qiH £ H 4x4 is a covariance matrix, then l’ q g = 1'^ H induces symmetries. As a consequence, there exists a more 
compact form of this matrix. 

Property 6: The covariance matrix of a centered quaternion random variable q takes the form: 


a 2 

7i^i 

7i/x. 

7]: Ml 

a 2 

71^3 1 

7im 2 

(7iM3 Ml )* 

a 2 




7i „ 3 


(71/xi A 


7i/x 3 

tW* 2 


where o 2 £ R, 7 i Ml £ H [1)M2)/ll3 ], 71 ^ £ and 7 i M3 £ H [1jMi>M2 ]. 

Proof: The term a 2 is real by definition. For 7 i Ml , using the Cayley-Dickson form q 
one gets that: 


(27) 


Zl + Z 2 P 2 with 2li ^2 £ C Ml , 


E [q( 5 Mi )*] =E[(a + /r 1 6 + /r 2 c + /r 3 d)(a-/r 1 6 + /r 2 c + /r 3 d)] 

= E [( 2:1 + Z2fJ-2)( z l + z 2P'2 )] 

= E [|zi| 2 - |z 2 | 2 + 2 Zl z 2 p 2 \ 

where 1 2 : 1 1,1 2 : 2 1 £ R. and Z 1 Z 2 P 2 * s a degenerate quaternion with only and parts, thus an element of H[ Mo ^ j. So, 
7 i Mi has no // | -part and is an element of By similar calculation, one can easily show that 7 i M 2 £ ] and 

7lM 3 ^ ® 

Property [6] shows how the knowledge of covariances of q with its three involutions gives the complete second order 
information about the quaternion random variable. It also demonstrates that the covariance matrix T q H is completely determined 
by one real number and three 3D degenerate quaternions, i.e. by a maximum of ten parameters. 

Expressions for 7i Ml , 7 i m , and 7 1/(; in terms of the previously mentioned Cayley-Dickson form of q are the following: 


<x 2 = E [|zi| 2 ] + E [|z 2 | 2 ] 

7 imi = ^[\zi\ 2 ] -E[\z 2 \ 2 ] +2E[ Zi z 2 \P2 

< (Z 0) 

7im 2 = E [ z i] + E [4] - 23f Ml (E [ 212 : 2 ]) p 3 
s 7 i„3 = E [4] - E [4] +( E M) P2 

In the sequel, we will sometimes make use of the following notation to avoid multiple indices: A = 71 , B = 7 lM2 and 

C = 7 i ; , for the covariances. These take their values according to: 


er 2 £ 


( ^ = 7i/x x £ e I[ 1 ,m 2 ,^ 3 ] ^ 9 ) 

B = 7 i M2 £ H [ ljAiij#i3 
C = 7in 3 £ ^[1,Mi, m 2 ] 

As previously mentioned, one can also study quaternion random variables with its real or complex vector representation, i.e. 
q : or q > Before inspecting symmetries for covariances of proper quaternion random variables, we give the expression of 
the covariance matrix using the complex vector notation as it will be of use in the interpretation of properness later on. First, 
recall using (|25]l and (p7| that with the quaternion vector representation, then: 
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V 2 

A 

B C 




r q ,n = Etqnqjj] 


A* 

a 2 

C Ml B^ 





B* 

(u*y 

a 2 






C* 


(A^ o 2 



Using the complex representation qc, the covariance matrix r q . :: takes the form: 




a 2 + R(A) 


+ Ql.Mil 

B [m 3 ]AG - <G/x 2 ]AG 

(^.MalA^) 

Tq,c = E[q c qc] = \ 

( S [l,Mil + C '[l.Mil) 

( 

/ 

7 2 + R(A) 

.. \ * 

A [m 2 ’M3] ^2 


“ B [m 3 ]AG - <G/x 2 ]AG 


( S [m 3 ]AG - C[h 2 ]Ag) 

[ A 

a 2 - K(A) 


( 5 ii.Md-qi.Mii)* 


A lZ 2 ,H 3 ]V 2 (' 

“ B [m : 

j] AG 


^Il.Md “ Ch.Mi] 

a 2 - R(A) 


where we used the relation F q c = 


GfcC 1 q,BM 


:C = 


with: 

1 
0 
0 


AG A*2 


AG AG 

0 0 


(30) 


(31) 


which relates complex and quaternion vector representations like: qc = Mitcqim- Recall that I’ q c takes values in C Ml which 
is isomorphic to the set of complex numbers, i.e. it is a complex valued matrix. 

In the next section, we investigate the consequences of the different levels of properness on the covariance matrix of 
quaternion random variables. 


B. Proper cases 

We now focus on different cases of properness and present the structure of the covariance matrix for some specific cases. 
Also, we point out the relation to previous work, especially the definitions from [3]. All the considered quaternion random 
variables from now are centered, i.e. IE//] = 0. There are actually four major types of properness which are of interest as they 
allow to express any proper case as we shall detail later. 

1) The {p 11 ii 2 )-proper case: This is the most “general” case. First, recall that (/x x , fi 2 )-properness stands for 7r ^ 2 (/x 1 , /r 2 ) 7r ^ 2 - 
properness. 

Definition 8: Given a quaternion basis {1, fi 1 , fi 2 . // 3 }, a quaternion Gaussian random variable q = Z\ + 22 AG with z\ . z 2 £ 
C M is called (/g, /r 2 )-proper iff: 


q = /g<7AG 


which involves the following symmetries: 


Tim, = 
/ l/x 2 


- 7V, 
71^2 

-Mi 


Mi 


Mi 


(32) 


(33) 


7im 3 — 7i„ 3 

The (/g, /x 2 )-properness has consequences on the covariance matrix elements, as stated in the following property. 

Property 7: Given a (p 1 , /x 2 )-proper quaternion random variable q, its covariance matrix, using the complex vector repre¬ 
sentation qc, takes the form: 


r q ,c = E[qcq, c ] = 


a 2 

a 


a 

a * 

a 2 

a* 

-aV 

-AG <5 

a 

a 2 

—a* 

a* 

AG<$ 

—a 

a 2 


(34) 
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where: 


<j = 


a = 


/V = 


E[N 2 ] =E[|z 2 | 2 ] 

E [zf] = E [ziz 2 ] = —E[^|]* 
E[^i z£] = -E[z^z 2 ] 


(35) 


a £ C Ml and (5 £ 


with (T 2 £ 

One can see from property [7] that a (/r,,, /i 2 )-proper quaternion random variable is completly described by four parameters. It 
can be also understood as a pair of improper complex random variables (z± and z 2 ) which are correlated and pseudo-correlatec 0 

Obviously, the properties axes and the way to write the Cayley-Dickson form of q ( i.e. the choice of Z\ and z 2 ) has 
consequences on the stmcture of the covariance matrix F q r. For example, (fi 2 ■ /t. 3 )-propcrness will lead to an other structure 
in T q c, but still only four parameters would completely describe the covariance. However, we do not give here the expressions 
of symmetries for 7 i m , and 7 i^ 3 , nor covariance matrices, for all possible pairs of axis (//. u) with // and u taking 

values n 1 , p 2 or /t, 3 , as they can be deduced by indices permutation in (33 i or computed directly using (30 1 . 

2) The (fi 1 , l)-proper case: This case is actually related to the definition given in [2], which makes use of left Clifford 
translation [18]. 

Definition 9: Given a quaternion basis {1, fi 2 , /r,}, a quaternion Gaussian random variable q = Z\ + z 2 p, 2 with z-\. z 2 £ 


C ;l is called (/j, ,. 1)-proper iff: 


d 

q = fitf 


(36) 


which involves the following symmetries: 


7 l,X, 

= 7ii*i 


7 i/x 2 

= -7iM 2 ' il 

(37) 

7im 3 

= - 7iM3 



The covariance matrix of the complex representation of a (/i, 1; 1)-proper quaternion random variable obeys the following 
symmetries. 

Property 8: Given a (/r ,, 1)-proper quaternion random variable q, its covariance matrix, using the complex vector represen¬ 
tation qc, takes the form: 


r q ,c = 


CT 2 0 CO 0 
0 CT 2 0 W* 


ur 

0 


0 ^ 2 

to 0 


(38) 


where: 


e[N 2 ] 
E [k 2 | 2 ] 

co = £[^ 1 ^ 2 ] 


C 2 = 


(39) 


with a 2 ,? 2 £ R. and to £ C Ml . 

From property [ 8 ] one can see that a ( p ,,. l)-proper quaternion consists in a pair of proper complex random variables with 
different variances (a 2 7 ^ ? 2 ), and which are correlated (E[ 2 i 2 2 ] 7 ^ 0 ) but not pseudo-correlated £[ 272 : 2 ] = 0 . 

3) The (1, ptfi)-proper case: This case is related to the definition in [3]. which is based on right Clifford translation [18], 
Definition 10: Given a quaternion basis {1, /x 1; /r 2 , /r 3 }, a quaternion Gaussian random variable q = Z\ + 2 2 /r 2 with 21 , 2 2 £ 
C Ml is called (1, ptfi) -proper iff: 

q = qp 2 (40) 

which involves the following: 

{ 7imi e 

7im 2 = 0 (41 > 

7i/x 3 = 0 


5 Given two complex centered random variables z\ and Z 2 , their cross-correlation (cross-covariance) is and their pseudo-cross-correlation (pseudo- 

cross-covariance) is £[ 212 : 2 ]. 
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The covariance matrix of the complex representation of a (1, p x ) -proper quaternion random variable obeys the following 
symmetries. 

Property 9: Given a (1, fi t )-proper quaternion random variable q, its covariance matrix, using the complex vector represen¬ 
tation qc, takes the form: 


r q ,c = 


(J 

0 

0 

w* 


0 0 

a 2 lo* 
oj < t 2 
0 0 


(42) 


where: 


(7 = 


E[N 2 ] 
E [M 2 ] 
u> = E[2i2 2 ] 


4 2 = 


(43) 


with cr 2 , ? 2 eR and 1 0 £ C /X| . 

From property [ 9 ] one can see that a (1, )-proper quaternion consists in a pair of proper complex random variables with 

different variances (cr 2 ^ <t 2 ), and which are correlated (E[zizJ] 0) but not pseudo-correlated £[ 212 : 2 ] = 0. 

4) The (p 3 , pf)-proper case: This last case of properness occurs when both axis are identical. 

Definition 11: Given a quaternion basis {1, p 1 , p 2 , p 3 }, a quaternion Gaussian random variable q = Z 1 + Z 2 P 2 with z\ , z 2 £ 

C, 


is called (p l7 p x )-proper iff: 


which involves the following: 


q = p ± qp ± 


(44) 


7iM! 

= 71(2^ 


7l(2 2 

= 7i(2 2 Ml 

(45) 

71(2.3 

= 7i(2 3 Ml 



The covariance matrix of the complex representation of a (p lt pf) -proper quaternion random variable obeys the following 
symmetries. 

Property 10: Given a (p l7 p ± )-proper quaternion random variable q, its covariance matrix, using the complex vector 
representation qc, takes the form: 


Tq,C = 


where: 


a 

a* 

0 

0 


cr = 


<T = 


0 

0 

? 2 

5* 


(46) 


a 

S 


E[M 2 ] 

E[|2 2 | 2 ] 

E[2 2 ] 

E[2 2 ] 


(47) 


with cr 2 , f 2 £ R. and with a,6 £ C 


Mi- 


From Property 10 one can see that a {p l7 pf )-proper quaternion consists in a pair of improper complex random variables 
with different variances (cr 2 c 2 ), and which are neither pseudo-correlated = 0) nor correlated (£[ 2122 ] = 0). 

Similar results could be deduced if we were to look at (p 2l /r 2 )-proper or (// 3 , /r 3 )-proper, and exact expression of covariance 
in those cases can be deduced from the presented results. This is also true for the three other properness cases: one can obtain 
similar results when replacing, in the described cases, the axes by p 1 , p 2 or p 3 . 

It is also interesting to note that, for second order study, there is no consequences in changing the sign of one of the axis in 
(p 1 , p 2 ) -proper definition. It means for example that (p l7 p 2 ) -proper and {p l ,-p 2 )- proper are equivalent. Other examples 
using any of the four levels of properness previously introduced can be derived. 
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C. Relation to previous work 

In [3], the three properness levels are defined using the vanishing of complimentary covariances. These covariances are in 
fact 7 i Ml , "n /,,, and 7 ^ as defined previously. The properness levels defined in [3] thus take the following expressions: 

• R-properness: 71 „ = 0 for one axis v> £ {/x x , /x 2 , /x 3 } 

. (C Ml -properness: 7 i M2 = 7 i M3 = 0 

K-propemesfj®] 7i Ml = 7i M2 = 7i M3 = 0 

O' 1 -properness from [3] is (1. p, { )-properness. The H-proper case in [3] is not an interesting 


Clearly, from Definition 


10 


one, as it simply consists in a quaternion random variable with uncorrelated components having the same variance. This case 
is actually a circular case in the Gaussian quaternion case (invariant by any rotation in 4D). The R-propcrness is actually a 
special case of (/r 1; /r 2 )-properness. As can be seen in [3], R-propcrness means that the quaternion random variable is a pair of 
improper complex random variables which are correlated but not pseudo-correlated. This is equivalent to being (1, -proper 
with additional constraint that Efzq z 2 \ = 0 (pseudo-covariance of Z\ and z 2 vanishes). 

In the approach introduced in [3], properness levels were strictly defined by the number of vanishing complimentary 
covariances ( 7 i Ml 71 1 , and 7-1 ;j;i ). It was done this way to mimic the complex case where an improper (resp. proper) complex 
variable has a non-vanishing (resp. vanishing) pseudo-covariance. However, in the complex case, the relation between properness 
and invariance by rotation of the probability density function is immediate. This is no longer true with the approach in [3]. 
With the (/i l5 /r 2 )-properness concept introduced here, the rotation invariances of the density can be directly related to the 
properness level, and consists in symmetries of the pseudo-covariances (see Definitions [ 8 ] [9 10 and m These invariances 
by (double) rotation of the quaternion random variable consist in a generalization of the work in [3] and [2], The properness 
definitions introduced in these references are thus special cases of p, 2 ) -properness. Actually, the definitions in [3] and [2] 
can not provide a complete picture of the possible symmetries for quaternion random variables, nor a way to explore circularity 
for those variables and its consequences on higher order moments and quaternion valued signal processing applications. 


V. Examples 

Examples of 4D distributions and constellations arise in applications such as optics (polarized signal processing) and have 
been studied in terms of optical devices random action (birefringence, etc.) in [27], [23]. Properness and improperness of 
quaternion random variables have also been exploited in filtering [8], detection [9], [10] and signal modeling [15]. Here, we 
simply illustrate the difference between the most currently used level of properness in quaternion signal processing [8], [6], 
[10], [15], [16], namely the C /x -properness, and the most general properness level defined in this article, namely the (/r 1; /im¬ 
proper. In order to simplify the notation, and without loss of generality, we present random variables expressed in the standard 
basis in H, i.e. { 1, i. j. k\. In figure [ 2 ] N = 5.10 4 realizations of a Gaussian (i.,/) -proper random variable are displayed. The 
4D plot is splitted into three pairs of 2D plots. Each pair of plots consists in viewing the 4D space as a pair of non-intersecting 
2D planes, namely: {l,i} and {;). A: }, ( l. j \ and ( i, k } and {1. fc} and [i, j }. The three plots are necessary to provide the 
complete 4D image of the possible correlations between the four components of the quaternion random variable. Figure [2] 
highlights the fact that for a (i,j )-proper Gaussian quaternion random variable, all the components are correlated. This is 
visible by the absence of any “proper” complex random variable (i.e. with rotational invariance) in any of the displayed 2D 
planes. In figure [ 3 } N = 5.10 4 realizations of a Gaussian (l,j)-proper quaternion random variable ar e display ed. This is the 


C J -proper case in [3] notation. One can see from figure [ 3 ] that in this case, and as explained in Section IV-B.3 the quaternion 


random variable consists in a pair of proper complex random variables (in planes {1, j } and { i, k\) which are correlated 
(improper plots in other planes). 

Property 11: Given a (l,j)-proper Gaussian quaternion random variable Q, then its probability density function is given 
by: 

1 ( 2o 2 \ q \ 2 -^(q*im j )\ 


fQ{q,q*,q 3 ,q ) = 


47r 2 det (T qiH ) 1/2 


exp -- 


-|7Gf 


(48) 


where det (r qj H) is the (/-determinant of T q> H as defined in [26], 

Finally, in order to illustrate the way (l,j)-proper has consequences on the expression of probability density functions, we 
give here the pdf expression in the quaternion vector representation <p-. One can see from Proposition that the pdf of a 
(l,j)-proper Gaussian quaternion random variable is only function of the modulus of q, i.e. 1 7 1, and of qq°, so that actually 
/q(9> <T l i W 3 >9~ fc ) = /q(<T 7~ j ) f°r a (l,j)-proper quaternion random variable. Similar reasoning can be transposed to other 
levels of (/x l 5 /r 2 )-properness with decreasing number of parameters as the symmetry constraints make the covariance matrix 
converge to a multiple of the identity matrix. 

Expression of the pdf of a (/i 1 ,/i 2 )-proper Gaussian random variable could also be derived using standard linear algebra 
calculus over H [28], and would consists in terms involving yu, 71 j and 71 Such expression could be usefull in properness 

6 Note that we use the notation H-properness here, as opposed to Q-properness in [3]. This choice is made to avoid confusion with the classicaly used 
notation for the set of rational numbers, i.e. Q, also used in this paper. 
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testing for example, and results similar to those presented in [6] could be derived. The properness test is in that case a LRT 
between hypothesis with different structured covariance matrices. Such study is left for future work. 

The use of (/x 1 , fi 2 )-properness in wiely-linear estimation algorithms should also have consequences and could lead to the 
development of new estimation algorithms (including the level of properness). Such algorithms could be designed for quaternion 
random vectors, for which the concepts presented in this paper can be straight forwardly extended. Finally, we emphsasiwe 
the fact that (/x x , fi 2 )-properness is a new tool to study pairs of complex valued variables/signals, and that it may be of 
practical use when considering quaternion valued random processes such as the stochastic version of the H-embedding signal 
(a quaternion-valued signal that can be associated to any non-stationary complex signal and allows a geometrical description 
of it) recently [29], 


VI. Conclusion 

A geometry based definition of properness allows to describe a wide range of different classes of quaternion random variables. 
Using the complete action of the rotation group SO{4) in R' 1 allows to identify levels of properness that where not known 
up to now for quaternion random variables. In particular, the geometric approach allows to identify “completely correlated” 
4D variables which have a certain level of properness and thus exhibit symmetries. The definition of (/r 1; /x 2 )-properness 
proposed in this work is a special case of circularity for quaternion random variables. The geometry (i.e. symmetries) of 
the pdf of the quaternion random variable involves a constrained structure on the second order moment (covariance matrix). 
(p, 1 , /r 2 )-properness is a more general concept than the previously introduced definitions, and actually incorporates these earlier 
definitions as special cases. Second order study of quaternion random variables can thus be performed in an exhaustive manner 
using our definition, and open the path to higher order statistics for such random variables. The use of p 2 )-properness in 

quaternion signal processing should allow to design new algorithms for estimation, filtering or detection that fully exploit the 
inner geometry of quaternion random signals. 
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Fig. 2. Realizations (N = 5.10 4 ) of a Gaussian (i, j )-proper quaternion random variable q. Four dimensional points are displayed in different pairs of 
non-intersecting 2D planes. Top: dl(q) vs ^Si(q) (left) and ^Sj(q) vs ^Sfc(q). (right) Middle: 9£(q) vs ^Sj(q) (left) and ^Si(q) vs Qfc(q) (right). Bottom: 9£(q) 
vs ^Sfc(q) (left) and $$i(q) vs ^Sj(q) (right). 
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Fig. 3. Realizations (N = 5.10 4 ) of a Gaussian (1, j) -proper, or C J -proper, quaternion random variable q. Four dimensional points are displayed in 
different pairs of non-intersecting 2D planes. Top: 3£(g) vs ^Si(q) (left) and ^Sj(q) vs (right) Middle: 3?(q) vs $Sj(q) (left) and ^Si(q) vs ^Sk(q) 

(right). Bottom: 9£(q) vs ^Sfc(q) (left) and $Si(q) vs ^Sj(q) (right). 

































































































